NASA Technical Memorandum 106120 


/ / \ / " • " . 

/ / v 

/7 

)i / ■' 

An Analytical Investigation of Transient Effects 
on Rewetting of Heated Thin Flat Plates 


J.A. Platt 

Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

Fundamentals of Heat Transfer in a Microgravity Environment, ASME 
1993 Winter Annual Meeting 

sponsored by the American Society of Mechanical Engineers 
New Orleans, Louisiana, November 28-December 3, 1993 


(NASA-TM-106120 ) AN ANALYTICAL 
INVESTIGATION OF TRANSIENT EFFECTS 
ON REWETTING OF HEATEO THIN FLAT 
PLATES (NASA) 11 P 


N9A-12968 


Unc 1 as 




G3/34 0179807 




AN ANALYTICAL INVESTIGATION OF TRANSIENT EFFECTS ON REWETTING 

OF HEATED THIN FLAT PLATES 


J.A. Platt 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

The rewetting of a hot surface is a problem of prime im- 
portance in the microgravity application of heat pipe tech- 
nology, where rewetting controls the time before operations 
can be re-established following depriming of a heat pipe. 
Rewetting is also important in the nuclear industry (in pre- 
dicting behavior during loss- of- coolant accidents), as well as 
in the chemical and petrochemical industries. 

Recently Chan and Zhang (1992) have presented a closed- 
form solution for the determination of the rewetting speed 
of a liquid film flowing over a finite (but long) hot plate 
subject to uniform heating. Unfortunately their physically 
unreasonable initial conditions preclude a meaningful analy- 
sis of start-up transient behavior. The current work pre- 
sents a new nondimensionalization and closed-form solution 
for an infinitely-long, uniformly -heated plate. Realistic 
initial conditions (step change in temperature across the 
wetting front) and boundary conditions (no spatial temper- 
ature gradients infinitely far from the wetting front) are 
employed. The effects of parametric variation on the result- 
ing simpler closed-form solution are presented and com- 
pared with the predictions of a “quasi-steady” model. The 
time to reach steady-state rewetting is found to be a strong 
function of the initial dry-region plate temperature. For 
heated plates it is found that in most cases the effect of the 
transient response terms cannot be neglected, even for large 
times. 


NOMENCLATURE 

b,c variables of integration 

c p specific heat 


h boiling heat transfer coefficient 

k thermal conductivity 

m function defined in equation (22) 

Q unif orm heat flux to the plate 

q dimensionless heat flux 

S plate thickness 

T temperature 

T 9 saturation temperature 

T 0 Leidenfrost temperature 

t time 

U rewetting velocity 

u dimensionless rewetting velocity 

X,Y length coordinates in Eulerian system 

X 7 length coordinate in Lagrangian system 

Greek symbols 

a,/3 functions used in temperature transformations 

(eqs. (11) and (18)) 

r) dimensionless length coordinate 

0 dimensionless temperature 

v transformed dimensionless temperature 

(eqs. (ll) and (18)) 

p density 

r dimensionless time 
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Superscript 


* 

reference quantity 

Subscripts 


d 

dry region 

w 

wet region 

wf 

wetting front 

1 

dry plate initial condition 


Y Liquid Y 1 

film — v 


— Wetting 
/ front 




Heat flux Q 

Figure 1 . — Geometry of the problem. 


INTRODUCTION 

Space Station Freedom (SSF), when built, will require 
more power than any previous space facility. The heat gen- 
erated as a result of this power must be transported to 
radiators where it will be discharged to space. One of the 
designs for the SSF thermal management system involves 
use of heat pipes in the radiators. Successful prediction of 
rewetting characteristics of heat pipes in a spacecraft envi- 
ronment requires an understanding of the physical mechan- 
isms that control the wetting of a hot surface. 

The problem of rewetting a surface which is significantly 
above the boiling temperature also has numerous industrial 
applications, most notably in the quenching of metals and 
in cooling nuclear fuel rod bundles during loss- of- coolant 
accidents in water reactors. Nuclear applications have been 
behind many of the previous analytical efforts at under- 
standing rewetting, such as those of Yamanouchi (1968); 
Thompson (1972); Duffey and Porthouse (1973); Coney 
(1974); Sun et al. (1974, 1975); Blair (1975); Tien and Yao 
(1975); and Dua and Tien (1976), which all found analytical 
expressions for the steady-state rewetting velocities by 
employing a variety of boiling models to various unheated 
surface geometries. 

Analysis of the unsteady start-up of rewetting and of the 
effect of uniform heating has been sparse. Bukur and Isbin 
(1972) performed a numerical analysis of start-up transients 
for a few cases corresponding to rewetting in nuclear loss-of- 
coolant accidents. They found that the effect of start-up 
transients for these cases was negligible. However, parame- 
tric ranges appropriate to heat pipe applications were not 
studied. 

A work has recently been published (Chan and Zhang, 
1992) which presents, among other things, a closed-form 
analytical solution of the unsteady rewetting of a heated 
thin flat plate. However, the initial conditions considered 
preclude any meaningful examination of start-up transients. 
Furthermore, no parametric analysis of plate heating was 
presented. 


In the present work a simpler closed-form solution for re- 
wetting velocity is derived for a thin heated flat plate of 
infinite length. Results are presented for a wide range of 
initial plate temperatures and heating rates, thus making 
the results applicable to microgravity heat pipe rewetting 
problems. The time-dependent start-up and heating effects 
of these solutions are discussed and full-solution heated- 
plate results are compared with those of a simplified quasi- 
steady solution. 


PROBLEM FORMULATION 

Consider the problem of rewetting, by a liquid film, an 
infinitely-long thin plate initially hotter than the Leiden- 
frost temperature (fig. 1). The plate is uniformly heated on 
its underside while the top dry region is considered adiaba- 
tic and the top wet region has a heat loss proportional to 
the difference between the plate temperature and liquid sat- 
uration temperature. The one- dimensional unsteady formu- 
lation of the problem, in dimensionless form, with 
appropriate boundary and initial condi tons, can be repre- 
sented as (see Appendix A): 


P° p 



dr 


_ u* u a O = k 

X* “ai") X* 2 Bf) 2 

+ S. -fi w (* < 0) 

(T 0 " T,)S S 


a* w 


dri 




= 0 


( 1 ) 

KM = 1 

KM = o 
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(v > 0) 



3«d 

3r 


U* u 9 ^d _ k ^ 9 A 
X* j X* 2 dr? 


(T 0 - T,)S 


d9 d 

a*? 


= o 




0 d (0,r) = 1 

* d M) = 0 1 


( 2 ) 


where r = , 

t* 





N| h h 


Thus the dimensionless systems 
and dry regions are: 


U* = 1 kh (6) 
pc p N S 

of equations for the wet 


3«w 

dr 


a 2 *,. a* w 

a, 2 a»? 


+ q - 


a* w 

w = o 
3* 


*w(M = 1 

* W M) = 0 


(»/ < 0) 


(7) 


(q = 0 is the location of the wetting front.) The reference 
time, length, and velocity are yet to be determined. In se- 
lecting dimensionless parameters the approach taken is to 
obtain unit-order dimensionless variables. This results in 
different parameters than those of previous investigators 
(such as Chan and Zhang (1992); Sun et al. (1974); and 
Duffey and Port house (1973)). One relationship between 
these reference values can be found by recognizing that the 
two terms on the left-hand side of the equation should be of 
the same order of magnitude, since they both arise from the 
unsteady term of the Eulerian formulation (eqs. (A4) 
and (AS)). Thus we have the relationship: 

U* as (3) 

t* 

Another relation is obtained by assuming that boiling is the 
dominant term on the right-hand side of the equation. (The 
assumption is amply justified for a problem involving 
sputtering.) 



Finally a reference length is chosen by balancing conduction 
with the unsteady and boiling terms: 




t* 



(5) 


d6 d 

"37 


a? e d 

~dn 2 



q 


39 d 




= 0 


*d(M = 1 


(r) > 0) 


( 8 ) 


9 d {lfi) = 


where q = 


Q 

S(T 0 - T,)‘ 


Here 0, the dimensionless temperature, is a function of 
position and time, and u, the dimensionless rewetting 
velocity, is a function of time. 0j and q are constants, with 
$ x greater than unity. 

In order to effect a closed-form solution to the problem, 
u will be treated as a constant for any given time. This 
approximation makes the equations linear and thus more 
tractable. Matching the rf = 0 temperature gradients for 
any given r will yield u for that time. 


SOLUTION FOR THE WET REGION 

First we solve the boundary condition as rj — > — «* to 
obtain an expression for 0 as a function of r. The relevant 
equation and initial condition are: 


Equations (3) to (5) can be combined to yield the appro- 
priate reference length, time, and velocity: 


— = q - U°) = o 

3r 


(9) 
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The solution is easily found to be: 


M- ->°) = *(! - e_r ) 


By reversing the transformation of equation (11) an expres- 
sion for the dimensionless temperature in the wet region is 
obtained: 


Note that for q > 1 complete dryout of the plate is expected 
for a finite time, ln(q/(q - 1)). These high heat flux 
conditions are excluded from further consideration. 

It is interesting to contrast the initial and boundary 
conditions of the present work with those of Chan and 
Zhang, who developed a closed-form solution for rewetting 
speed on a heated plate of finite (but very large) length. 
Chan and Zhang’s boundary condition to from the wetting 
front was 0(-i? Ll ,r) = 0 (where i? Ll » 1), which models the 
end of the plate in contact with a liquid reservoir at 
saturation t emp erature and results in spatial temperature 
gradients far from the wetting front when q # 0. 

Chan and Zhang’s initial condition for both the wet and 
dry regions was 0(*?,O) = $ v This is a physically unreason- 
able condition for the wet region and as a result the 
rewetting velocity could not be determined for small times. 

By employing a transformation to get the problem in the 
form of a conduction equation: 

». » [#. - ,d - 

u u 2 

where a w = _ and /3 W = 1 + — 

2 4 

the following system of equations is obtained: 

= 0 

(12) 

uJO,r) = 1 - q(l - e->^ r 

v w( - ".») = 0 


This system is very close to the elementary conduction 
problem of a semi-infinite body, initially at constant tem- 
perature, that is subject to a step change in boundary tem- 
perature at time 0. By employing Duhamel’s theorem 
(Myers (1987)), we can use the solution to that problem to 
obtain a solution for the present problem as: 




ZJL f 1-1(1 


4(r-c) 


dc 


( 13 ) 


= <i(l - e T ) ~ 


tie 






x j (i - 3)T^ + 


c=0 


(’ ~ 


x exp 


4(r - c) 


dc (14) 


This can be rewritten as: 

fi w = q(l - e- r ) + (1 - q)e- (a - 4 ^ ) ’’ 
+ q e-e- (a - 4 ^ n_) ’’ 




4- 


f J 

£7 b =° 


(1 - q) exp - 


i? 2 b 2 


+ qe' 


. r f (?* + ')* 

exp — — 

l b 2 


„ 2 b 2 

4t 


b 2 4f 
db (15) 


where b 


f 




Substituting the expressions for a w and j3 w , the wet-side 
temperature gradient at t) = 0 is 


M w 

dr. 


- !Vu 2 + 4 - u)(l - q) + i^u 2 + 8 - u)qe r 


4=0 


+ 2(l^q) j exp 


/xr~ 


b=o 


_ (u 2 + 4)r 


4b 


db 


4- 


! exp 

Vxr b =° 


(u 2 + 8 )t 


4b 


db 


(16) 
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SOLUTION FOR THE DRY REGION 

The solution for the dry region follows a s imil ar procedure 
to that for the wetted region. In this case the boundary 
condition far from the wetting front is: 

*d(“> T ) = 6 i + q T ( 17 ) 


Upon applying the transformation: 

v d = (*d - V - 6 1 )e ad ” +ffdT 


(18) 


where = _ and 


2 u 

= Qj = — 

d 4 


the following system of equations is obtained: 
^d _ 

= o 


v di.°’ T ) = (1 - *1 - q* , )e* <r 

*'d(°". T ) = 0 


(19) 


By application of Duhamel’s theorem the following results 
are achieved: 


9 d = K + q r ) - ( 9 i + qr - i)e 2 “ a ” 

+ (0 1 + q t — i)e adV ^ 


/ 7TT 


dO d 

'*T 


1 

x f exp 

Pd T »j 2 b 2 

bio 

b 2 ?) 

2u(0 L + qr - 

1) 


db 


■• 1=0 




1 

c 

_fr) 

exp 


b =0 

, 4b 2 J 


( 20 ) 


db (21) 


Matching expressions in equations (16) and (21) for the 
temperature gradient at t} — 0 yields an expression that, 
given $ 19 q, and r, can be solved iteratively for u: 


2u(^ + qr - 1) 



+ 4 



- q) 



+ 8 




x [(*i + q r “ l)m(u 2 ) - (1 -q)m(u 2 + 4) 

- qe _r m(u 2 + 8)] = 0 ( 22 ) 


where m(x) = J exp 


b = 0 


XT 


4b A 


db 


This expression is significantly simpler than the only other 
comparable solution in the literature, equation (52) of Chan 
and Zhang (1992). 

For q = 0 and r — ► «, equation (22) reduces to a solution 
equivalent to the one first derived by Yamanouchi (1968): 


u 


N 


4 

2{e 1 - i ) 2 - i 


(23) 


RESULTS AND DISCUSSION 

The dependence of velocity on initial plate temperature, 
for the case of unheated plates (q = 0), is shown in figure 2. 
As expected, the rewetting velocity is reduced as the initial 
plate temperature is increased. 

An unusual characteristic of the present solution is the 
initial infinite rewetting velocities for $ x + 2. This is caused 
by the temperature step-function initial condition, which 
results in infinite wetting-front temperature gradients in 
both the wet and dry regions as r — ► 0. The sign of these 



Dimensionless time, t 

Figure 2. — Velocity versus time for unheated (q = 0) conditions. 
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infinite velocities depends on the relative temperature 
differences between the wet and dry initial temperatures (0 
and 0 X , respectively) and the wetting-front (Leidenfrost) 
temperature, 0=1. 

Note that this step function initial condition, while pro- 
ducing infinite positive or negative initial velocities, does 
not produce an infinite initial movement of the wetting 
front. This can be shown by examining ur for small r. 
Unless ur — ► 0 as r — ► 0 the physically impossible result 
occurs of a nonzero movement of the wetting front, either 
finite or infinite, for vanishingly small time. Such a physi- 
cally impossible result does not occur for the present 
solution (fig. 3). 

A significant result of the present work is that the time to 
reach steady-state velocity increases sharply as the tem- 
perature of the plate increases (fig. 4, solid line). For 0 1 = 2 
about lOt* is required to reach steady state while for = 5 
about 100t* is required. 

The implications of this result for nuclear applications are 
minimal because the high boiling heat transfer coefficients 
for such cases result in small reference times. Using the val- 
ues of density, specific heat, and plate thickness given by 
Bukur and Isbin (1972) for the experiment of Yamanouchi 
(1968), and using a heat flux of 10 6 kcal/m 2 hr °C 
(Yamanouchi, 1968), t* = 0.035 sec. (The numerical solu- 
tions of Bukur and Isbin (1972) confirm that the effects of 
initial conditions are minimal for nuclear applications.) 

The time to reach steady state can be significant for other 
applications, for which the boiling heat transfer coefficient 
is several orders of magnitude lower. For the experiment of 
Grimley et al. (1988), using values derived by Chan and 
Zhang (1992), t* = 7.4 sec and the time to reach steady 
state is approximately 150 sec. Thus for nonnuclear situa- 
tions such as heat pipe rewetting, the amount of time to 
reach steady state could be significant. 

Results for heated plates (figs. 5 to 7) establish that, 
given enough time, for all cases velocities will first become 
negative and then will approach zero. (The dry region plate 
temperatures for which these transitions occur are so large 
that in most real situations they would be physically unrea- 
sonable.) The fact that u — + 0 as r — ♦ «> can be accounted 
for simply by the effect of increased dry region temperature; 
as the temperature of the plate gets very large rewetting 
occurs very slowly. 

An explanation of the negative velocities is more difficult 
and can be aided by a comparison of the results of the 
complete equation (eq. (22)) with those of a “quasi- steady” 
solution based on just the first two terms of equation (22): 

F /2 

u = 4(1 - q) 2 (24) 

[(20, + 2qr - 1 - q) 2 - (1 - q) 2 

The velocity in this quasi-steady solution is unconditionally 
positive (the same would also be true for a solution based 


— Time to reach 99% of steady-state velocity 

— 2 + q t (plate temp, when 8-j = 2) 



Figure 3. — Time to reach steady-state velocity (solid line) and 
plate temperature tor 0-j = 2 heated cases (broken lines). 



Dimensionless time, t 
Figure 4. — ut versus t for small t. 


on the first three terms of eq. (36)), indicating that the 
source of the negative velocities is in the last three terms, 
system. 

A comparison of the full solution and the quasi- steady 
solution for B x = 2 (fig. 8) confirms the import ance of the 
transient terms. Agreement between the two solutions 
occurs only for q = 0.01. The reason for the lack of agree- 
ment in the other cases is evident from an examination of 
figure 4, where the broken lines map the plate temperatures 
for B x — 2 and various levels of heating. Only the q = 0.01 
curve ever reaches what can be termed a “quasi-steady” 
state, where the velocity can be accurately approximated 
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Figure 5.— Effect of heat flux on velocity, 8-| = 2. 



Figure 7.— Effect of heat flux on velocity, O-j = ^ .5. 


by an infinite time solution at the same plate temperature; 
and even in that case the solution becomes unsteady again 
for larger values of r. Therefore it is concluded that, for 
q > 0, the effect of the transient terms is significant in 
almost all cases. The negative velocities are evidently the 
result of localized heating in the wetting front region, as 
opposed to the effect of heating on the dry plate 
temperature far from the wetting front. 


CONCLUSIONS 

An analytical investigation has been undertaken of the 
transient rewetting of an infinitely-long, heated flat plate 



Figure 6.— Effect of heat flux on velocity, = 2.5. 


Full solution 



Figure 8. — Comparison of full closed-form solution with "quasi- 
steady" solution, 0-j = 2. 


initially above the Leidenfrost temperature. An approxi- 
mate closed-form solution for the rewetting velocity has 
been obtained by applying some simplifying assumptions to 
the equations for heat conduction in the plate. 

With a step function initial temperature condition, the 
initial rewetting velocity is infinite (either positive or neg- 
ative), although this result is explainable and is not physi- 
cally unrealistic. The time required to reach steady-state 
wetting velocity is found to increase sharply with increasing 
plate temperature. 

The effect of heating the plate is to decrease rewetting 
velocity over time, with dryout (negative rewetting 
velocity) eventually occurring in all cases (although the 
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plate temperatures for which this dryout would occur would 
generally be too large to be achieved in an actual system). 
“Quasi-steady” states (heated conditions for which transient 
terms do not affect the results) are found to occur for only 
a narrow range of time and heat flux. 


APPENDIX A: DERIVATION OF EQUATIONS 

For the problem of rewetting an infinitely-long plate ini- 
tially hotter than the Leidenfrost temperature (fig. l) the 
two-dimensional) unsteady conduction equation for tem- 
perature in the plate is: 


nx^t) = T 0 





(X > X wf ) 


3T 

ax 


i(^) 


= o 


(A5) 


KX^t) = T 0 


' C P 


3T 

at 



3Y 2 , 


(Al) 


The plate is uniformly heated on its underside while the top 
dry region is considered adiabatic and the top wet region 
has a heat loss proportional to the difference between the 
plate temperature and liquid saturation temperature. Thus 
for a thin plate the conduction term in the Y-direction can 
be simplified as: 


k 


a*T 

BY 2 


Q 

S 

Q 

s 



T.) 


( wet > (A2) 


(dry) 


Far from the wetting front, conditions will be unaffected by 
movement of the wetting front and temperatures will be 
independent of X, leading to boundary conditions of: 


3T = Q (A3) 

3X U-.t) 


Upon transformation to a Lagrangian coordinate system 
(X',t) with its origin at the wetting front equations (A4) 
and (A5) become: 


dT TT 3T 

' Cp lar U ax; 


= ki!Z. + S-h(T - TJ (X'<0) 
3X' 2 s 


3T 

ax 1 


= o 


(A6) 


(--.*) 

T(0,t) = T 0 




[aT_u3T) =k i!l + S (X' > 0) 

l at ax'J ax- 2 s 


3T 


= 0 


(A7) 


aX 'l(-,t) 
T(0,t) = T 0 


Therefore we have two one- dimensional, unsteady conduc- 
tion problems to solve (one each for the wet and dry 
regions). An additional boundary condition for each of the 
problems is the temperature at the wetting front 
(X = X wf ), which is assumed to be the Leidenfrost temper- 
ature. The systems of equations (without initial conditions, 
which are not yet specified) are: 


PS 


aT 

~dt 



h ( T _ T j 

s 


(X < X wf ) 


aT 

ax 


l(--kO 


= o 


(A4) 


where U is the rewetting velocity, which in general is a 
function of time. 

Initial conditions are sought which satisfy time- 
independent, U = 0 forms of equations (A6) and (A7). For 
the dry region an isothermal initial condition is the only 
possibility: 


T(X' > 0,0) = T x (A8) 


where T 1 is a temperature greater than the Leidenfrost 
temperature. For the wet region, although another solution 
is possible, an isothermal initial condition, with the plate at 
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the saturation temperature of the liquid, is adopted for 

simplicity: 

T(X' < 0,0) = T, (A9) 

The initial conditions therefore are nothing more than a 

step function across the wetting front. 
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variation on the resulting simpler closed-form solution are presented and compared with the predictions of a 
“quasi-steady” model. The time to reach steady-state rewetting is found to be a strong function of the initial dry-region 
plate temperature. For heated plates it is found that in most cases the effect of the transient response terms cannot be 
neglected, even for large times. 
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